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Thermodynamic Magnon Recoil for Domain Wall Motion 
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We predict a thermodynamic magnon recoil effect for domain wall motions in the presence of temperature 
gradients. All current thermodynamic theories assert that a magnetic domain wall must move toward the hotter 
side, based on equilibrium thermodynamic arguments. Microscopic calculations on the other hand show that a 
domain wall can move either along or against the direction of heat currents, depending on how strong the heat 
currents are reflected by the domain wall. We have resolved the inconsistency between these two approaches 
by augmenting the theory in the presence of thermal gradients by incorporating in the free energy of domain 
walls by a heat current term present in non-equilibrium steady states. The condition to observe a domain wall 
propagation toward the colder regime is derived analytically and can be tested by future experiments. 

PACS numbers: 75.30.Ds, 75.60.Ch, 85.75.-d 


Spincaloritronics is the subfield of spintronics which ex¬ 
plores spin-dependent phenomena coupled to thermal gradi¬ 
ents 0], A very important question within this field is how 
a magnetic domain wall (DW) can move under a tempera¬ 
ture gradient. This question has attracted much attention ow¬ 
ing to its applicability in magnetic insulators for poten¬ 
tial applications in logic devices J3l and data storage tech¬ 
nology iH]. The conventional approach using static mag¬ 
netic fields II] is well established with high DW velocities 
IToliTil . but does not allow for the synchronous motion of 
multiple domain walls. Synchronous current-induced domain 
wall motion due t o sp in-transfer torques 131415tl and/or spin- 
orbit torques 116142011 formed an alternative way to efficiently 
manipulate the magnetization configuration, but the required 
high-current densities cause problems, such as Joule heating 
owing to Ohmic losses. Heat itself has been proposed as an 
efficient control parameter to overcome the problems during 
the emergence of spincaloritronics 12111 . Fully understanding 
and predicting new controlled ways to move the domain walls 
by magnonic heat currents is paramount to exploiting fully all 
their future device possibilities. 

There is at present an theoretical incomplete understand¬ 
ing of temperature-gradient-driven DW motion. There are two 
types of theories, i.e., a macroscopically thermodynamic the¬ 
ory M and a microscopically magnonic one I 22 I - I 24 I 1 . The 
theories contradict each other in certain regimes. In previous 
thermodynamic theories M, a magnetic domain wall at fi¬ 
nite temperature T is treated as a thermodynamic object with 
free energy, F = U - TS , where U is its internal energy and 
S is its entropy. The free energy of the DW can also be ex¬ 
pressed as the difference between a system with a DW minus 
that of the same system without the wall A F - A U - TAS. 
Thermodynamic calculations M show that, far below the 
Curie temperature, the entropy AS ( T ) increases and the free 
energy A F (T) decreases with the temperature. This leads to 
a conclusion that the DW must move towards regions with 
higher temperatures due to the entropic force, with the prop¬ 
agation velocity proportional to the temperature gradient WT 
12 - 3 ]. This tendency has been observed in experiments fl]- 

On the other hand, there are microscopic angular momen¬ 


tum transfer 122] and linear momentum transfer 12311 theo¬ 
ries for the magnon-driven DW motion. The DW moves 
along the opposite direction of current flows if the former 
mechanism dominates S [zl] . while along the same direc¬ 


tion when the latter one is more important 12314250 , i.e., there 
are strong spin-wave reflections by the wall |23,j2f|. The 
proposed mechanism of magnonic linear momentum transfer 
has been_confirmed in various systems including ferromag- 
nets 123, 24l l, antiferromagnets lf27[|28l l, and spin textures with 


Dzyaloshiskii-Moriya interaction ]29 

In this Letter, we show that the inconsistency between these 
two types of theories can be resolved by augmenting the free 
energy by a term from the heat current which always flows in 
a non-equilibrium steady state in the presence of a tempera¬ 
ture gradient. The heat current gets modulations by the DW 
with momentum-conserving back-scatterings, ft then leads to 
a force that pushes a DW toward the colder region. We pre¬ 
dict a new thermodynamic magnon recoil effect for the do¬ 
main wall motion in temperature gradients. Under conditions 
of a strong backscattering, a high magnon thermal conduc¬ 
tivity, and a slow magnon group velocity, this magnon recoil 
effect surpasses the previously identified entropic force §3- 
Such a regime can be achieved in yttrium iron garnet (Y1G) 
and other ferromagnetic insulators, as we show below. 
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FIG. 1: (Color online) Heat currents in a ID magnetic wire in the 
presence of a temperature gradient VT without (upper panel) and 
with (lower panel) a domain wall with width <5 fl , denoted by j’ and j^, 
respectively. The presence of the domain wall causes the reflection 
of heat currents with a probability R. 
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We consider a one dimensional (ID) magnetic wire con¬ 
necting two thermal reservoirs with temperatures 7'/, and T\ 
( T ], > Tj) as shown in Fig. 1. The temperature gradient V7' 
then drives magnon heat currents jj and jj? in a uniformly 
magnetized wire (as shown in the upper panel in Fig. 1) and 
in a wire with a domain wall (shown in the lower panel in Fig. 
1), respectively, at non-equilibrium steady states. The pres¬ 
ence of the domain wall causes finite magnon reflections 0, 
thereby reducing the heat current, i.e., jj? < jj. They are equal 
only when the potential generated by the domain wall is re¬ 
flectionless |23 23]. A heat current jj 11 can modify the rate of 
the change of the entropy |30| 31 j 


dS 1 ’ 11 

dt 




( 1 ) 


and this entropy change modifies the free energy. We there¬ 
fore write the free energy of the wall as 


A F = A F e + AF n 


( 2 ) 


Here A F e is the equilibrium part as treated in previous thermo¬ 
dynamic theories ffl , which may originate from the modi¬ 
fied magnon density of states due to the wall 00El, while 
AF ne is the non-equilibrium part proportional to the heat cur¬ 
rent. Its change rate is then 


d(AF ne ) 

dt 



where we introduced a parameter R(T) denoting the 
temperature-dependent reflection probability of magnon heat 
currents by the wall. The magnon heat current in the pres¬ 
ence of the domain wall is thus modified, i.e., jj 1 = (1 - R) j*. 
The heat current is connected with the temperature gradient 
by jj = -kVT, with a positive-definite magnon thermal con¬ 
ductivity k. We then obtain 


d (AF„ e ) 
dt 



(VT) 2 dx > 0, 


( 4 ) 


for finite reflections of heat currents. The temperature depen¬ 
dence of the free energy A F (T) is crucial to drive the DW 
propagation. It has been shown that the equilibrium part A F e 
decreases with an increasing temperature T (10], which leads 
to a conclusion that the DW must move to the hotter region to 
reduce the free energy due to the entropic force. However, 
the present system under the temperature gradient V7\ strictly 
speaking, is not in an equilibrium state. Moreover, the non¬ 
equilibrium part A F ne monotonically increases with time [for 
R > 0 in Eq. Q I, thus the arguments based on either max¬ 
imization of entropy or minimization of free energy are not 
always valid [2H4!] since transport processes matter a lot in 
non-equilibrium thermodynamics. In the following, we study 
the magnon transport, in particular its backscatterings by a do¬ 
main wall, and predict a thermodynamic magnon recoil effect 


in the presence of a temperature gradient, competing with the 
entropic force. 

The thermal properties of magnons crucially depend on 
their dispersion relations and lifetime. To calculate the 
magnon thermal conductivity k, we consider the heat current 
carried by the magnon flow due to the temperature gradient 
V7' in the absence of domain walls (as the upper panel in Fig. 
1 ), jj. = L Yjk 6nkhu)\ g (k), where L is the wire length, k is 
the magnon wave-vector, 6 tik = tik - hi_ is the magnon num¬ 
ber in excess of equilibrium value hk = 1/ | e M*)/(*nD _ ij 
being the Bose-Einstein distribution with Boltzman constant 
kg , hu (k) is the magnon energy, and v„ (k) = da>/d k is the 
magnon group velocity. Using the Boltzmann approach we 
can write a first-order expression for the excess magnon num¬ 
ber in the steady state and in the relaxation time approxima¬ 
tion, 5nu - -Tk (dFik/OT) v,,-V7\ where t/ c is the magnon relax¬ 
ation time. One thus obtains the magnon thermal conductivity 

/r=—V Tkhu> (dfik/dT) Vgdu>, (5) 

ln 


by using the one dimensional magnon density of states (DOS) 
p(k ) = L/2n. Here N is the number of energy bands and 
wjjun(max) j s lowest (highest) frequency of each band n. 

The presence of domain wall may lead to a strong spin- 
wave reflection, and thus a reduction of magnon heat cur¬ 
rents jjj = (1 - R) jj.. The reported mean free path of ther¬ 
mal magnons in insulating ferromagnets, e.g., YIG, usually is 
~ 1-100 fi m |0, 34j], which is much larger than the domain 
wall width 6 d ~ 10 - 100 nm, the scattering of spin waves 
by the wall can thus be treated as a ballistic process, thereby 
conserving the total momentum. We then can derive the re¬ 
flection probability R of magnon heat_currents by the wall via 
the Landauer-Buttiker formula 1 0 35] 


R(T) = 


Zn=i£f F(co, T) \r(k)\ 2 du> 

XLjX F(0J ' T ^ d0J 


( 6 ) 


where r(k) is the /.--dependent reflection coefficient of 
magnons by the wall and F ( a> , T) — tia> ( dhk/dT). Here we do 
not consider the modification of the magnon DOS due to the 
wall 12611 . which is relevant to reflectionless magnons treated 
in equilibrium thermodynamic theories [0, 32 ] but causes only 
negligible effects to our results here. In the momentum- 
conserving scattering process between spin waves and the 
domain wall, the change rate of the linear momentum of a 
DW is dpuw/dt = 2cj)M s ly 123ll which must be compensated 
by that from magnons (with wave vector k ) dp mdgnons /dt = 
(Srik) v g |r (k )\ 2 tik. Here cf> is the tilted angle of the DW plane, 
M s is the saturation magnetization, and y is the gyromagnetic 
ratio. Spin-wave reflections thus lead to a precession of the 
domain wall plane with the angular velocity 


h = ( 5n k ) v g |r (k)| 2 hk. 


( 7 ) 


The equivalent magnetic field responsible for the above pre¬ 
cession velocity is then h* = <j>k/y^ giving rise to an effective 
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field along the wire axis after a summation of all states 
H ne =L-'^h, 


vr 

AjiM, 


N 


;v r>i 

Zf, 

0 ), 


Tk (ditk/dT) Vo |r (k)\ 2 fikdoj , (8) 


which is the effective field or force acting on the wall due to 
the thermodynamic magnon recoil in a temperature gradient. 

Equation ® is quite a general formula that can be used to 
calculate the effective field under any magnon dispersion re¬ 
lations and relaxation mechanisms. Experiment data in YIG, 
for instance, show an acoustic branch with frequency that rises 
from nearly zero at the Brillouin zone center to a value at the 
zone boundary that varies from 6 to 9.5 THz. These values 
correspond to temperatures of approximately 300 and 500 K. 
Since the lowest optical branch lies above the zone-boundary 
value, the calculation of the thermal properties up to room 
temperature can be done considering only the acoustic branch, 
i.e., N — At low wave numbers the dispersion relation 
can be approximately by a quadratic form w = w™ n + Jk 2 , 
where w™ n is the acoustic band gap depending on materi¬ 
als parameters, such as the magnetic anisotropy, dipole-dipole 
coupling, Dzyaloshiskii-Moriya interaction, etc., and J is the 
exchange constant. The cut-off frequency w™ ax thus is a)" 13 ' = 
w™ n + Jk ^ with k m the maximum wave vector depending on 
the magnon propagation direction. The magnon group veloc¬ 
ity is then v g - 2-Jj (w - , 


u>. 


It has been shown that the 


quadratic dispersion agrees very well with the actual disper¬ 
sion up to a wave vector k = 0.6 k m in YIG [36|]. Under the 
above conditions, we obtain 


H ne = -- 


kVT \r\ 2 
2nM s Vo ’ 


(9) 


with the average reflection probability 


i -|2 


r u T 

JuT n 


TkF (to,T)\r(k)\ 2 dcj/ I t^F (to,T) du>. 


r"i 


( 10 ) 


and the average group velocity 

TkVgF (a>,T)da>/ I TkF (w, T) da>. (11) 

w min mn 

We then obtain the DW velocity along the direction of heat 
currents due to their recoil effect, below Walker breakdown 

v„en = ^H„, = Ji£2 (i2) 

a v g a l7iM s 

Schlickeiser et al. 10] derived an effective mag¬ 
netic field based on equilibrium thermodynamics, H e = 
VT (JqIT c ) / (<i D>n eq M s aJ with nearest-neighbor exchange en¬ 
ergy Jo, Curie temperature T c , equilibrium local magnetiza¬ 
tion m eq , and lattice constant a, which is believed to exceed the 


magnonic spin transfer torque proportional to 1 — |r| 2 fHH. 
So, the final DW propagation direction depends on the com¬ 
petition between H ne and H e . The condition to observe a DW 
propagation toward the colder regime is therefore 


2 n 


v g 6 D m eq a 


(Jo/Tc), 


(13) 


which requires a good heat conduction in magnetic domains (a 
large k without the domain wall), a strong magnon backscat- 
tering (a large \r | 2 ), a slow magnon group velocity (a small Vg), 
and a broad domain wall (a large 6d)- The forces induced by 
non-equilibrium thermal fluctuations under temperature gra¬ 
dients cause a Brownian motion |37j, 38] of the domain wall 
and could be another reason to push its propagation toward the 
colder region 09]. However it is still an open question how 
valid the classical fluctuation-dissipation theorem for equilib¬ 
rium states is S, particularly when it is applied to nonequi¬ 
librium steady states in the presence of temperature gradients 
141 -44]. The Brownian motion effect is however negligible in 
the presence of strong magnon backscatterings. 

In order to evaluate the parameters in criterion (IT3l> . we 
now can make either of two plausible assumptions about 
the behavior of Model I: If one considers that the re¬ 
laxation time Tk is independent of both the wave number 
and the temperature (the simple average-lifetime model) 14511 
and takes t> = f k, we obtain |r[ 2 = R the same as the re¬ 
flection probability of magnon heat currents [Eq. <[6]) when 
N = 1], Vg = v g F(oj, T) doj/ F (w, T) doj, and 

k oc f k- Model II: If we consider the Gilbert damping 
but neglect higher order processes such as magnon-magnon 
interactions, the relaxation time is then t> = 1 / (2aa>) 

with Gilbert damping constant a 140]. We thus have 
lif - tfTm/d T )\r(k)\ 2d aj/ tfT m k /dT)daj, v g - 

J(x) j ° 

max max 

(dnk/dT)vodu>/ f (dhk/dT)du>, and k oc 1/a. 

6 J(jJ\ 

It has been shown that both dipole-dipole 12311 and 
Dzyaloshiskii-Moriya (29] interactions can result in strong 
magnon reflections in the presence of a domain wall in fer- 
romagnets. A precessing domain wall in antiferromagnet can 
also lead to significant magnon reflections 128, 47]]■ A com¬ 
mon feature of the reflection probability function |r(k(«))| 2 
is the sharp transition from 1 at lower frequencies to 0 at 
higher frequencies II23Ll24i 14711 . satisfying ansatz |r(k(«))| 2 = 
|w(- (u> - a> c ) /Aw)| 2 , with the transition frequency a> c and 
the spectrum width Aw. Function |w| 2 reduces to 1 for a> <s 
a> c , and 0 for o> » oj c . The form of the function w depends 
on material parameters such as the domain wall width |[26f| , 
the Dzyaloshiskii-Moriya interaction strength |29i], etc., and 
scattering details such as the incident angle of magnons [23]. 
However, for a very narrow spectrum (Aw « 1) which is of¬ 
ten the case [23l 14711 . it can be approximately described by the 


Heaviside step function, i.e., |w| 2 « s (- (w - w c )) with 

( 0, x < 0 
s (x) = i i, x — 0 . 

1, x>0 


( 14 ) 




























4 


It indicates that magnons are completely reflected by the wall 
when their frequencies are lower than a> c , while they es¬ 
sentially pass through the domain wall without any reflec¬ 
tion above u> c . By denoting A = /k B , B = ki<jj™ m /k B , 

x a = A/T , and x c = C/T with C = ha> c /k B , Eqs. (0, (fTOt . and 
(flTT i can be calculated analytically and yield 


Tk 


Jk\T 3 


rB/T 

Ja/t 


y/x - x a x 2 e* 


v„ — 2 


Irl 2 =R = 


p __ 

Jk B T Ja/t (e*-if 

h f (A/T) - f (B/T) ’ 
B * T s(—(x — x c )) , e , dx 


'a/t ( e x - 1)" 

■B/T VFFTA 


-dx 


dx. 


Ia 


<a/t 


(e*-ir 


f (A/T)-f (B/T) 


with |48fl 


f(x) = -2x In (1 - e~ x ) + 



(15) 

(16) 
(17) 


(18) 


in Model I (t> = Tk), and 


1 Jk B T C 

2 na y n J A 

- r 

Jk B T Ja 


Jk 3 B T r B,T y/x - x a xe* 


Vg - 2 


A/T 
B/T yjx-x a xe 
'A/T {e x -l) 


(e x - 1) 
-;-dx 


H g(A/T) - g (B/T)’ 


C s ' (~ (x ~ x c )) 


Irl 2 


JA/T 


(e x —l)' 


t dx 


g (A/T) - g (B/T) 


with |4g] 


dx, 


(19) 


( 20 ) 


( 21 ) 


g(x) = - In (1-0 + -^-, (22) 

e x - 1 

in Model II (j k = 1/(2 aw)). 

The following three cases are of potential interest, (i) For 
a low transition frequency (oj c < w™”), there is no reflection, 
and Eqs. ( fTTI) and (l2ll yield f| 2 = 0. (ii) For an intermediate 
transition frequency (w™ n < a> c < w™ ax ), Eqs. (fT7] > and (EH 
reduce to |r| 2 = [/ (A/T) - / (C/T)] / [f (A/T) - f (B/T)\ 
and [g (A/T) - g (C/T)] / [g (A/T) - g (B/T)], respectively, 
(iii) For a high transition frequency (a> c > w™ ax ), all magnons 
are reflected by the wall. Thus, Eq. < fl7T> and (EH reduce to 
|r| 2 — 1- 

Figure 2(a) shows the temperature dependence of the 
magnon thermal conductivity k in both Models I and II. It 
increases as the elevated temperature in two cases. In Fig. 
2(b) we calculate the average group velocity v g . It mono- 
tonically increases with the temperature and saturates at high 
temperatures. The temperature dependence of |r| 2 for differ¬ 
ent transition temperatures C are shown in Fig. 2(c), with a 
monotonically decreasing manner. It is because higher tem¬ 
peratures make more magnons populate higher energy levels. 



T(K) 


T(K) 


FIG. 2: (Color online) (a) Temperature dependence of the magnon 


heat conductivity k normalized by T k 



in Model I, and 


by (2 Tia) 1 yjjC B A/h 


in Model II, respectively. 


(b) Average group 


velocity v g as a function of temperature T , in unit of 2 y/Jk B A/Ti 
in two models, (c) Temperature dependence of average reflection 
probability \r\ 2 for different transition temperatures C. (d) Parame¬ 
ter k |r| 2 /v g as a function of temperature for different C, with unit of 
f k k 2 B A/ (2nTi) in Model I and of k B / (4 na) in Model II, respectively. 
In the calculations, we use A = 100 K and B = 500 K. All tempera¬ 
tures are below T c . 


which leads to a smaller magnon reflection subsequently. We 
also observe that a higher transition temperature leads to a 
stronger magnon reflection. Figure 2(d) demonstrates a mono¬ 
tonically increasing dependence on the temperature of param¬ 
eter k\ r\ 2 /v g in both models. As shown in Eq. ( fXH with 
quadratic magnon dispersion relations, parameter k f| 2 /v g , 
which is independent of exchange constant J, is crucial to de¬ 
termine if the domain wall can move toward the colder region. 
According to our calculations, this condition should be easily 
satisfied at elevated temperatures (T > 400 K) with a small 
magnon damping ( f k > 1 ns or a < 10 4 ) for a broad domain 
wall (6 d k. 100 nm) in a weak ferromagnet (Jo/a < 3 x 10~ 12 
J/m) under any temperature gradient that can overcome the 
pinning force produced by defects or impurities. Other relax¬ 
ation models considering three- and four-magnon scattering 
processes are expected not to modify our conclusions signifi¬ 
cantly. 

To summarize, we predict a thermodynamic magnon recoil 
effect for domain wall motion under temperature gradients. 
We correct the previous thermodynamic theories by includ¬ 
ing a heat current term for entropy and/or free-energy gen¬ 
erations, which is always presents in non-equilibrium steady 
states in the presence of a temperature gradient. The heat cur¬ 
rent gets modulations by the DW with momentum-conserving 
backscatterings. It then leads to a recoil force on the wall, 
which competes with the previously identified entropic force. 
Our theory thereby closes the inconsistency between macro- 
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scopic and microscopic theories for the domain wall motion, 
and and we propose experiments to test it. We also expect the 
similar thermodynamic magnon recoil effect to play an impor¬ 
tant role in other magnetic structures, e.g., magnetic vortices, 
bubbles, or Skyrmions, and other materials like antiferromag- 
nets or multiferroics. 

After the completion of this work, we became aware of one 
recent report [42j] on a DW thermophoresis in antiferromag- 
nets using the classical fluctuation-dissipation relation without 
considering any magnon backscattering. Our results should 
also be applicable to their work. 

We thank Se Kwon Kim and Gerrit Bauer for useful dis¬ 
cussions. This work is supported by DFG Priority Pro¬ 
gramme 1538 “Spin-Caloric Transport” (P.Y. and J.S.) and 
European Union Seventh Framework Programme “Spinlcur” 
under Grant No. FP7-People-2012-ITN-316657 (Y.C.). 
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